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In the paper we investigate the theory of quantum optical systems. As an application we integrate 
and describe the quantum optical systems which are generically related to the classical orthogonal 
polynomials. The family of coherent states related to these systems is constructed and described. 
Some applications are also presented. 
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I. INTRODUCTION 

Quantum optics affords a big amount of very inter- 
esting physical phenomena having important application 
at the same time. Our aim is to formulate the theory of 
these phenomena and elucidate their connection with the 
theory of orthogonal polynomials. This allows us to use 
the last one for the rigorous integration of some nonlinear 
quantum optical models describing the interaction of the 
finite number of modes of electromagnetic field with non- 
linear medium. Let us mention, that in quantum optical 
literature, see e.g. §, §, 0, (l5| the solutions 

of the models of this type are usually approached by ap- 
proximative or semiclassical methods. 

In Section 2 we deduce from natural and not restrictive 
assumptions, the general form, see ( p.35| ), of the Hamil- 
tonian H/ describing the interaction of the finite number 
of modes of electric field with the mater. Later on, in 
Section 3, we investigate the quantum reduction method 
which allows us to describe the quantum optical systems 
by the use of the theory of orthogonal polynomials (see 
0, d], Q). We also show that the reduced syste m ar e 



related to some quantum algebras, see relations < ^.27\ ,. 
These algebras were investigated in ||2^, where their re- 
lations to the theory of special and g-special functions 
was shown. 

Having spectral measure of the interaction Hamilto- 
nian, which is for example the case if the model under 
consideration is related to the classical orthogonal poly- 
nomials or their g-deformation, see we can in- 
troduce spectral coherent states. They are direct gener- 
alizations of Glauber coherent states (corresponding to 
the Hermite polynomials) and squeezed states. In Sec- 
tion 4, we show that spectral coherent states admit the 
holomorphic representation for the Hamiltonians under 
consideration which supplements the spectral and Fock 
representations. This simplifies remarkably the calcula- 
tion of many important physical characteristics of the 
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described system. 

The spectral coherent states should have some physical 
meaning which needs the deeper understanding. In any 
case, they give the link of quantum optical systems with 
complex analytic and symplectic geometry. This opens 
the application of coherent states method, investigated 
in [gO|, 1 21 to the problems of the theory of quantum 
optics. 

In Section 5. we give complete solution of the quan- 
tum systems, see (pT|), related to the classical orthogonal 
polynomials. 

Finally, in Section 6, we present th e phy sical interpre- 
tation of the Hamiltonian given by (2.37) as the para- 
metric modulator, which includes as special cases such 
quantum optical systems as nondegenerate parametric 
amplifier and the frequency up-converter, see [ p7| . 

At the end we express our conviction that the pro- 
posed method will be helpful in better understanding of 
quantum optical problems. 



II. QUANTUM ELECTROMAGNETIC FIELD 
IN NONLINEAR MEDIUM 



Nonlinear optics deals with phenomena that occur as a 
consequence of the modification of the optical properties 
of a material system in the presence of light. Practically, 
only laser light is sufhciently intensive to produce the 
measurable effects. By an optical nonlinearity we mean 
that the dipole moment per unit volume, or polariza- 
tion P, of a material system depends in nonlinear way, 
upon the strength of the applied electromagnetic field. 
As many authors js), ||2^ we assume that P depends 
only on the electric part E of the electromagnetic field 
{E,B) i.e. P = P[E]. We assume moreover that this 
dependence is a functional one. Thus, in most general 
case we can write 

P[E] (t, f) = £0 ^ / dti (ffi ...dtN d^TN (2.1) 

\n) {t,r,ti,ri, ...,tN,rN) (^E{ti,ri) , E {In ,rN)j , 
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where the vector valued A'^-hnear map T'(Ar), is called in 
optical literature iV-th response tensor of the medium . 

The time-invariance principle, which says that the dy- 
namical properties of the system are assumed to be un- 
changed by a translation of the time origin, leads to 



T{N}{t,r,ti,ri, 
-■ R(N){f,ti 



■ ,tN,rN) 



,ijv-i,rjv).(2.2) 



The interaction of the electromagnetic field {E, B) with 
a nonlinear medium characterized by polarization P can 
be described by the source-free Maxwell equations 



ot 

V- (eo^ + P[^]) =0, 
V • B = . 



(2.3) 



Therefore the divergence of the Poynting vector -^^Ex B 
takes the form 

J-V.(i?xS) = -If^S^ 
Mo ot \2^Q 



2 



^ dP ^ 



where the quantity 



Uq 



1 

2^0 ' 



-B' 



£0 



E' 



(2.4) 



(2.5) 



is the energy density of the free electromagnetic field. 
Analogously, we define the interaction energy density ui 
by the equation 

The energy density 

u{t,f} ■.= uo{t,r)+ui{t,r) (2.7) 
determines the Hamiltonian H of our system 

H = Ha + Hi, (2.8) 



where 

i/o = / Uo{t,r}cfr 



(2.9) 



is the Hamiltonian of the free electromagnetic field and 
the Hamiltonian 



Hi = y" ui{t,r)(Pr 



(2.10) 



describes the interaction of electric field E with the 
medium under consideration. 

In order to obtain an explicit formula for ui let us 
consider the electromagnetic field potential A 



efe,A^A(fc)e^("'=*-'-"^^ 



(2.11) 



expressed in terms of Fourier modes, where the index 
A € {1,2} labels the polarization of the field, which is 
described by the pair of unit vectors Ck^i and or- 
thogonal to the wave vector k (we choose the Coulomb 
gauge V • A = 0) . Here we do not specify the form of the 
dispersion relation, so we assume that ujk is any function 
of |fc|. In this gauge we have 

d 



B = y X A. (2.12) 

Let us introduce the following simplifying notation 

-.a _(-iuJkek.x forcr=l 

^ - for(7 = -l ^^■^'^> 

^J(fc):=(^^(?) for<J=l 

' \ Al{k) for a = -1 ^ ' 

We have now 

E{t,r)^Y.j Al{k) e-("'=*-^"-^) d'k (2.15) 

and, therefore, (EJ) becomes 



oo „ 

P[E]{t,r)^^ ^ d^ki . . .(fk^ X{N) {r,ai,ki,uJki, ■ ■ ■ ,crN,kN,(^k,. 

AT n _ \ _ \ " 

N 

it (Tr-UJkr- 



N—Q cri,Ai ctn^^T' 

N 

'Allik,)...AZik^), 

where XiN) ci, fci, Wfe^ , . . . , trjv, kN,i^kN^ , is the 7V-th susceptibility tensor defined by 



X(N)(r,(Ji,ki,ujki,--,<JN,kN,(^kN^^ R{N) {r,Ti,fi,...,TN,rN) e ==i 



(2.16) 



dn d^fi...dTn d^rN-i^.U) 
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Inserting (2.16) into (2.6) we find up to additive constant that 

CO „ 

N 

J2 <^r^kr 



Ui( 



TV— cro,Ao cri,Ai (Tj^jXj^ 



r=l 



N 
s=0 



(2.18) 



In the quantization procedure the classical quantities are moreover, replaced by the normally ordered products 
A^(fc) in ( ^Hl) are replaced by the operators of corresponding operators, i.e. by : 

In orde r to obtain the Hamilt onian it is e noug h to in- 
(2.19) sort (|2.19| ) into ( ^.18|) and then (|2.18|) into (|2.10| ). With 
the Hamiltonian of the free electromagnetic field 



^k,\ 



^k,\ 
*fc,A 



for (7 = 1 
for (T = — 1 



which satisfy the commutation relations of a free quan- 
tum field: 



^A,/c J^A'.fc' 



^AA'<5(fc-fc')2(l-^^')- (2.20) 



The products A^°(fco) . . . A^^(fcAr) in ( |2.16| ) and ( |2.18| ) we obtain 



Ho = XI / '^'^ ^fe,Aafc,A 



(2.21) 



CX) „ 

H = Ho+X X ... X j(fkQ...(fk_ 

AT — n _ \ _ \ 



N ■ '^koM-- kN^N ■ ^ 



it 2 a^Uks 



AT— OtTo,Ao criv,A 



where 



AT 

' - fci,c^fc,,..., (TAT, ^Ar,c^fcJe-'^«*^«-^^"dV. (2.23) 



(- 



s=0 



Using the commutation relations ( 2.20 ) one can prove that 



iHot „(T iHot _ pirrujkt a 



k,X 



Hence the Hamiltonian (2.8) becomes 
H = Ho + e-'"°' Hi e'^°\ 



(2.24) 
(2.25) 



where due to ( p.22| ) 

oo „ 

Hj = X X . . . X d^ko... d^kis 



I 



does not depend on time, and therefore the solution of is given by 

the Schrodinger equation 1^^^^^ ^ ^^^IIot^-^u,t j^^q^^^ 

'|l*<'» The operator 



(2.28) 
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'iHot 



(2.29) 



is the free electromagnetic field evolution operator. The 
operator 



Vi{t) 



-iHit 



(2.30) 



is the evolution operator of the system in the interaction 
picture. 

For the real models in quantum optics one assumes 
that the system under consideration contains a finite 
number of modes of electric field (see , , iH , 12^ 



). This means that the label (fc, A) in ( |2.26| ) and (2.21 
takes a finite number of values 



(fc,A) = je {0,1,... ,M} 



(2.31) 



and the integrals are reduced to finite sums over j: 



Ho = E 



N=0 (To ,jo <tn Jn 



(2.32) 
(2.33) 



In this case the Hamiltonian ( 2.33 ) can be transformed 
into the form which is more useful for our aims. It is 
defined by the exchange of the normal ordering of the 
annihilation and creation operators into the one which 
we will call boson-number ordering in the sequel. 

In order to define the boson-number ordering let 
us introduce the following notation for creation and an- 
nihilation operators 




for Z = 1,2, . 
for 1^0 
for / = -1,- 



(2.34) 



A product of m annihilation and n creation operators be- 
ing in the same mode, is said to be boson- number ordered 
if it is of the form P(a*a)a™~", where P is a polynomial. 

Changing the normal ordering in each term of 
the Hamiltonian ( |2.33| ) to the boson-number ordering 

al^/ where P is a polynomial 



To 



of M -|- 1 variables and collecting the terms with the same 



factor a. 



H, 



we obtain 



... JmSZ 



9la,... ,1m (^0^0' ■ 



)a[,V. 



(2.35) 



where gig^...j,^j are functions of (M -I- 1)— variables de- 
pendent on 0(7v)- The Hamiltonian H/ is a symmetric 
operator if 

[gio,... :1m (^0^0' ■ • ■ ' ^M^Af )]* 

9-i„:... :-Im i^Q^o -[o,-- - , a*Ma.M - hi) ■ (2.36) 
In the next sections we restrict our considerations to the 
Hamiltonians of the form 



H/ — h (aQaQ, . . . , a^^a^^^) 

+.9 (^0^0' ■■ ■ ' ^M^Af) ^0 • ■ ■ ^ 



(2.37) 



Such form of the Hamiltonian is strictly related to the 
theory of orthogonal polynomials. The physical inter- 
pretation of this Hamiltonian is given in the Paragraph 
6.A. 



III. REDUCTION OF THE HAMILTONIAN 

In this section, we briefly describe the decomposition 
of the Hilbert space H spanned by elements of the or- 
thonormal Fock basis 



Bf 



\no, 



, Um) 



.nM\ 



{air 



no,... ,nM € NU{0} 



(3.1) 



into invariant subspaces of the operators Hq and H/. 
The method of this decomposition is presented in details 
in j22j. In such a way we obtain the reduction of the 
Hamiltonian H. 

The invariant subspaces of Hj are obtained in two 
steps. The first step is related to some family of inte- 
grals of motion; the second one is related to a family of 
pseudo-vacuum vectors. 

Let us start with a few definitions: 



A/ * ^ \ Iq ' 

■= 9 l^-QaQ, . . . , a^a^ j ag ... a 



M 



(3.2) 



and 



M 

E"u a*a^ , 

3=0 



z = 0,l,... ,M, 



(3.3) 
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where a = (ay) is a real (M + 1) x {M + l)-matrix 
satisfying the conditions 



det a ^ 0, 



M 



(3.4) 
(3.5) 



The invertibihty of the matrix a ahows one to express 
the boson- number operators SL*a^ by A^, which gives 



H 



M 
3=0 



(3.6) 



with real constants jj determined by the matrix a. In 
particular we have 



M 
i=Q 

Additionally 

Hj = Hd(Ao,Ai, 



(3.7) 



(3.8) 



with Hd uniquely determined by the function h and the 
matrix a. Using the canonical commutation relations for 
creation and annihilation operators one obtain 



AA* -g(Ao,Ai,... ,Am), 
A*A = g(Ao-l,Ai,... ,Am), 
[Ao,A]=-A, [Ao,A*]=A*, 
[A,,A]=0, j = l,...,M, 

[A„A,]=0, i,j = 0,...,M 



(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 



with the nonnegative function Q uniquely determined by 
g and a. 

Direct calculations gives 

[A,,Ho] = [A,,Hj] = j = l,2,...,M. (3.14) 

which means that operators Ai,A2,... ,Ajv/ are inte- 
grals of motion. 

In order to reduce Hq and H/ to the common 
eigenspace of integrals of motion let us notice that the 
operators A* A, AA*, Aq, . . . , Am are diagonal in the 
Fock basis Bp- This, in particular, means that each vec- 
tor |no, tim) € is the eigenvector of the operators 
Aj, j = 0, . . . , M , with eigenvalues given by: 



M 



1=0 



(3.15) 



Moreover, the operators Aq, . . . , Am form a system of 
commuting independent observables. In such a way we 
can use the sequences of eigenvalues (Aq, Ai, . . . , Am) as a 
new parametrization {|Ao, Ai, . . . , Am)} of the Fock basis 
elements. So we obtain 



Since [Aq, A] = —A then, from (3.10) and (3.9) we have 
A |Ao, Ai, . . . , Xm) 

= (Ao - 1, Ai, Am) |Ao - 1, Ai, Am),(3.17) 

A* |Ao, Ai, . . . , \m) 



= V5 (Ao, Ai, Am) |Ao + 1, Ai, Am) . (3.18) 

It is clear that the subspace 'H\^...Xjy, of the Fock 
space H spanned by the eigenvectors |Ao,Ai,... , Am) 
with fixed Ai,... , Am is Ho and H/-i nva riant and 
dimH\j^...\nj = oo if and only if all Ij in ( |3.2[ ) are non- 
negative. The problem of integration of the system ( 2.25| ) 
is reduced to integration of the system described by the 
reduced Hamiltonian 



M 



H 



0,red 



7oAo + ^7jAj, 



H/^red Hd (Ao, Ai, . . . , Am) + A 
and therefore (up to additive constant) 



7oAo ■ 



ed e 



170 Ao* 



(3.19) 
(3.20) 

(3.21) 



Now we go to the next step of the reduction. In order 
to make it let us define the pseudo- vacuum vector as such 
vector |Ao, Ai, . . . , Am) from the Fock basis in TYai,... .Am 
which is annihilated by the operator A, i.e. 



A|Ao,Ai,... ,Am) =0 
or equivalcntly 

g(Ao-l,Ai,... ,Am) =0. 



(3.22) 



(3.23) 

In ||2^ it was shown that the set {Xo,i}h=i lA p : 
A|Ao,Ai,... ,Xm) = 0} of the solutions of ( 3.22 ) is 
nonempty if in the definition ( |3.2| ) any Z^, i = 
0, 1, . . . , M, is greater then zero. 

Now, if for simplicity, we introduce the notation 



b{n) 
N 

then 



= |Ao,i + n, Ai, . 

= y^^(Ao,i +n 
= Ao — Ao,; , 



, Am), 



l,Ai, 



, Am), 



(3.24) 



N|n) — n\n), 
A\n) = 6(n)|n-l), 
A*|n) = b{n + l)\n + l). 

Thus we obtain that the space 

J- :— span{ \n), n = 0, 1, . . . } 



(3.25) 



(3.26) 



is the irreducible representation space for the algebra 
Ared generated by the operators N, A and A*, which 
satisfy the relations: 



Aj |Ao, Ai, Xm) — Xj |Ao, Ai, Am) , J = 0, 



[N,A] = -A, [N,A*] = A*, 
A*A - b\N), 
(3.16) AA* = 



(3.27) 
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These algebras were investigated in [^. The question 
when the dimension of T is finite or infinite was discussed 
in detail in |22|. Here we assume that dimT = oo. After 



restriction to JF, the Hamiltonians ( 3.f 9 ) ( 3.2C ) belongs 
to Ared and take the form (up to additive constant) 



Ho,red — 7oN, 



A*, 



(3.28) 
(3.29) 



where h{N) := Hd(N + Aq,;, Ai, . . . , Am)- Thus the op- 
erators Ho. red , H/^red and Consequently Hred belong to 
the algebra Ared- In the Fock basis { \n), n = 0, 1, . . . } 
the operator H/^red assumes the three diagonal (Jacobi) 
form: 

= h{n)\n) + b{n)\n - 1) + b{n + l)|7i + 1), (3.30) 



whereas Hp, red is diagonal 

Ho^redl") = 70"|")• 



(3.31) 



From now on we restrict our consideration to the space 
J^. In particular we restrict all operators discussed above 
to and omit the index red for simplicity. 

The evolution of the system given by ( ^.2^ ) now takes 
the form 



and therefore for any operator F we have 



(3.32) 



{m\Fm) 

= (V'(0)|e*"^*e*^''^*Fe-'^«^*e-'«^*|i/;(0)). (3.33) 

For a special, but interesting case this formula simplifies. 
Namely, 



m)\fmm) 

= (^(0)|e*H^7(N)e-^"^*|V'(0)), 



(3.34) 



= /(e-*^''*)(V(0)| e^"^V(A) e-*H^*|^(0)), (3.35) 

imiUm) = 7o(V'(0)|e^«^*Ne-^«^*|V^(0)) 

+ (^(0)|H,^(0)), (3.36) 

where / is an analytic function. 

The next section is devoted to the detailed study of the 
operator H/ and 1-parameter group e^*^^* generated 
by it. 



IV. SPECTRAL AND COHERENT STATES 
REPRESENTATIONS 



The operators H/ of the type ( p.2£ ) are very well 
known in the theory of orthogonal polynomials 



. They are symmetric in J- and, by ( |3.30 ) , have a dense 



domain which consists of finite linear combinations of el- 
ements of the Fock basis. The deficiency indices of H/ 

oc 

are (0,0) or (1,1). One can prove that if ^ = oo, 

n— 1 

then the operator H/ has deficiency indices (0, 0), which 
is equivalent to its essential selfadjointness. 

From now on we will assume that the deficiency indices 
of H/ are (0,0). Hence H/ admits a unique selfadjoint 
extension, which will be denoted by the same symbol. 
Moreover, H/ has simple spectrum. This fact allows us to 
identify the Fock space with the Hilbert space of square 
integrable functions L^(M, da) of real variable w € R. 
The measure da is determined by the spectral measure 
dE of the hamiltonian H/ and is defined by the formula 



da{uj) 



'\dE{uj)0). 



(4.1) 



Additionally one can prove, that polynomials {uj"'}'^^q 
form a linearly dense subset in L^(M, da). After the 
Gram-Schmidt orthonormalization of the basis {w"}^o 
we obtain an orthonormal set {Pnj^o in L^(]R, da) 
called the orthonormal polynomial system. Notice that 
deg Pn = n. 

The unitary isomorphism U : J- — > I? (M, da) of 
Hilbert spaces is given by 



U|V^) :=^(n|^)P„ 



(4.2) 



Tl = 



According to the spectral theorem and (4.2) one has 

(Uo/(H,)oU-i)V(cc>) = /(co)^(co) (4.3) 



for ij ) g £^ (R, dcr) and any mea surab le function /. By 
(4.2) and (O) the expression ( 3.30 ) converts into the 
three-term recurrence formula 



wP„(w) = h(n)P^{uj) + fe(n)P„_i(w) -I- h(n -|- l)P„+i(w) 

(4.4) 

for the system of orthonormal polynomials {PnjJ^o- '^^ 
by the spectral theorem in the notion of this orthonormal 
system, we have: 

(m|/(H/)n>- f f{Lu)P^{uj)Pr,{Lo)da{Lo). (4.5) 



In particular, for /(H/) = -pjHj we obtain the moments 



(4.6) 



fik of the measure ( |4.1| ) : 

Mfe / u;' daiu) = ^{0\ HJ|0). 



Similarly, for /(H/) = we obtain the absolute 

Ml 

moments \fik\ of (4.1): 



I \u;\''daiu;)^^{0\ |H,|'= |0). 



(4.7) 
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oo 
k=0 



For the case under consideration the moments {fik} 
determine da in the unique way 

From (4.3) and (4.5) one obtains that the evolution 
operator e^^^*, t e M, in the Hilbert space L^(R, dcr) is 
given by 



(U o e"*"^* o U-i) ^/.(w) 



V(^) 



(4.8) 



and its mean value in the vacuum is realized by the char- 
acteristic function 



e-*"*dcr(cj) = 



1 



" |0) (4.9) 



of the measure da, compare with |10| . 

After these preliminary remarks we will show that 
apart from realizations of the Hamiltonian H/ in the Fock 
space T and in the Hilbert space (M, da) it is useful 
and natural to consider its realization in some Hilbert 
space which consists of square integrable holomorphic 
functions defined on an open subset of complex plane. 
To do it let us first prove the following 

Lemma IV. 1 Let us assume that absolute moments |/i|„ 
are finite for all n G NU{0} and they satisfy the condition 



lim - 

n — ^cxj 



eR 



< +00. 



(4.10) 



Then there exists a maximal strip in <C, which is open, 
connected and invariant under the one-parameter group 
of translations 



Ttz :^ z + t, t (z 



(4.11) 



such that the characteristic function a{t) can be holomor- 
phically extended to it. 

The maximality of the strip means that a{t) cannot be 
extended to a larger set with the same properties. 

Proof : We prove firstly that characteristic function a 
is analytic on the strip \Imz\ < R. One has 



dt" ^ ^ ^ ' 



'u;''da{u) 



(4.12) 



for n = 0, 1, . . . . In order t o pro ve ( 4.12 ) one proceeds 
by induction. The equality (4.12) is valid for n = 0. Let 
us assume that it is true for n, then 



p-CT(t) = lim 



^ — i{t-\-h)Lj 



lim 



h 
- 1 



-a;" da{uj) 



to 



da{uj) 



Since 



ihuj 



n+l 



(4.13) 



(4.14) 



and |/x|„ is finite for n = , 1, . . . , we were able to use 
Lebesgue theorem in ( 4.13 ). For S R we have the 
estimate 



fe=0 

a{t + h) - 



k\ dt^ 



E 



i-ihuf 



k=0 

E 

fe=0 



kl 



' da{uj) 



fc! 



where for the last inequality we used 



S fe! 

fc=0 



< 



By the Cauchy criterion and Stirling formula 



i! = V2^ n" e"" e®(") 



(4.15) 



(4.16) 



(4.17) 



where Q{n) < , the series ^ ^fI/^I" convergent 



for \h\ < R. This and ( 4.15 ) imply that Taylor expansion 

°° jn 

^(^ + ^)-E^^^W (4.18) 

fc=0 

is convergent for \z\ < R. We have proved the analyticity 
of a on the strip \Im,z\ < R. So there exist a nonempty, 
maximal strip {z e C : 2r < Imz < 2s}, such that 
the characteristic function a(t) can be holomorphically 
extended to it. We have shown, moreover that this strip 
contains the real axis i.e. — oo < 2r < < 2s < +oo. 
QED 

Let us consider a "half of the strip {z e C : 2r < 
Im z < 2s} i.e. 



Yi :— {z (z C : r < Imz < s}. 



(4.19) 



As a consequence of the Lemma IV.l we can formulate 
the following 

Prop osition IV. 1 Under the assumptions of Lemma 
IV. \ the map 



K -.^3 z< — > e-'^ ■ e L^(R,da) 



(4.20) 



is holomorphic and its image K (S) is linearly dense in 
L^{R,da). 

Proof : 

In order to see that the function e~*^' belongs to 
i^(K, da) let us notice that 



da{uj) = j e-*(^-^")" da{uj) 

= a{z - z) < +00 (4.21) 
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for (z — z) S {z S C : 2r < Im z < 2s}. Thus in the basis 
{Pn}n=o we have 



(4.22) 



n=0 



where the coefficients functions (t„ are holomorphic ex- 
tensions of 



(4.23) 



onto the whole strip S. Thus the map K is complex 
analytic map of the strip E into Hilbert space (R, da) . 
In order to show that K{T,) is linearly dense in 
(K, da) let us notice that the monomials 

^" = *"^^WMU=o, (4.24) 

where n = 0, 1, . . . , belong to the linear closure of if (S). 
Since, they form linearly dense subset of L^(M, da) and 
the same property is shared by K^S). QED 



Combining (Q with (|^) we obtain a holomorphic 
map 



n=0 



dn{z)\n)eT (4.25) 



of S into Fock space JF. Following ||2^, we shall 
call K : T, ^ the coherent states map related to the 
quantum system described by the Hamiltonian H/. The 
states |z), where z G E, will be called spectral coher- 
ent states. The coherent states map has nice physical 
properties and, as we will show later, it is useful for the 
calculations of physical characteristics of the system. 
By the formulae 

ft :— i-——(\oga{z — z))dz/\dz 
azaz 

1 

= --^{\oga{2iy))dx^dy, (4.26) 

where z = a; + iy, we will define the symplectic form 
on E. Using the mean value function 

(z|H/|2) _ 1 d 



loga{2iy) 



{z\z) 2 dy 

of the Hamiltonian in spectral coherent states 
we define the classical Hamiltonian system 



(4.27) 



z) ,z e E, 



(4.28) 

on the symplectic manifold (E,ri). The Hamiltonian 
flow, tangent to the vector field X^Hj)^ is given by (4.11). 



Let us denote by CP{J-) the complex projective Hilbert 
space modelled on the Fock space T. Let ilps denote 
Fubini-Study (1, 1)— form on CP(J-'),(for the definition of 
rips consult for example [0). The form Qps is closed 
and nonsingular. So, (CP(J-'), fi^s) can be considered 
as a symplectic manifold which can be interpreted as 
the quantum phase space of the system described by the 
Hamiltonian H/. 



Proposition IV. 2 The projectivization JC : S — > 

CP{J-), )C{z) := C\z) of the coherent states map {{.2b) 
is the holomorphic symplectic map, i.e. 

IC*nFs^^ (4.29) 
and the diagram 



E ^ CP(J^) 
E ^ CP{T) 



(4.30) 



is commutative for any t G 
Proof 



The equality ( [4.29| ) can be checked by direct calculation. 
The c ommutativity of the diagram (4.30) follows from 
( 4.22 ) by the use of the formulae for quantum (L8) and 
classical (4.11) evolution of the system. QED 



Recapitulating: we see that the coherent states map 
maps symplectically the classical phase space E of the 
system (E,f2, (H/)^) into the quantum phase space 
CP(JF) of the system (CP(JF), 51^5, H/). It is equivari- 
ant with respect to the classical and the quantum flows. 
The mean value function (H/)^ of the quantum Hamil- 
tonian H/ give the classical Hamiltonian of the system. 
So, the above picture is analogous to the one related to 
the harmonic oscillator (see pj|). For the general the- 
ory of quantization and description of physical systems 
in terms of the coherent states map see |^ . The model 
of the physical system considered here gives an important 
and interesting example illustrating the theory which was 
developed in pi[ . 

Let us define spectral annihilation operator a by 
the condition 



a z 



z\z), 



(4.31) 



which means that a has the spectral coherent states |z) 
as eigenvectors with eigenvalues z € E. It is defined on 
the dense linear domain, spanned by spectral coherent 
states. The representation in P^(M, da) is given by 



(U o a o U^V) {^) = i-ri'ioj)- 

■ ' dw 



(4.32) 



The domain D (U o a o U^^^ 
of all polynomials. 

According to Proposition 



is given as a vector space 



IV. 1 



the spectral coherent 



states form a linearly dense subset in T . Hence one can 
define anti-linear monomorphism U 



T3\i,)^ U|V^) := (^.|if(.)) e 0(E) 



(4.33) 



of the Fock space T into vector space ©(E) of holomor- 
phic functions on E. In such a way we obtain the third 
realization of Hilbert space of states, this time as the 
space of holomorphic functions U(JF) C C(E) with the 
scalar product defined by 



(cI>|vI/)^(U(|0)),U(|V))) := m). 
where ^I^ = Ul-^), $ = U|0)- 



(4.34) 
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Proposition IV. 3 Let the measure 
d^j,{z, z) = ^{y)dxdy, 



(4.35) 



on S, (z — X + iy) be such that the weight function ^ 
satisfies 



da ^ 
dcu 



M / dy ^,{y)e'y'^ = I 



(4.36) 



for ui £ supp da. 

Then the scalar product (4-34) can be expressed by the 
integral 



(4.37) 



(^10) = / ^{z)^z)d^i{z,z). 

Moreover, the kernel function 
{z\v) — a{v — z) 



(4.38) 



is a r eprod ucing kernel function with respect to the mea- 
sure { 4--3t ), i.e. 



-^(v) ^ I a{v - z)'^{z) dfi{z, z) 



for any * G U(J^). 



(4.39) 



Proof : 

In order to prove that ( 4j7| ) and (4.39) are valid for 
dfi{z, z) given by (4.35 - 4.361 ) let us observe that (i^(z, z) 
has the form (4.35) since the kernel a{- — z) is invariant 
with respect to the one-parameter group of translation 
(4.11). Hence we have 



a{y-z)a{z-w)dii{z,z) = da{uj) da{uj') / dx ^(y)dj/ e-*"'^+™'^'e^"("'-'^)e^("+'^') 



da{u:) / da{uj')5{uj - uj') / ^^/^(j/) e«('^+'^')e"™'^+™'^' 



da[Lo) f ^{cu + t)S{t) dT dyii{y) e!/(2-+-)e-»(— -Vg™" 
J did 

da{Lo)^{co) f dyii{y) ^-'^-^^^ . 



(4.40) 



If p satisfies (4.36) then ( 4.40| ) takes the form 

'a[y ~ z)a{z — w) diJ.{z, z) = a{v — ui), (4-41) 



which the reproducing property. QED 

In the sequel, let us assume that U(JF) — i^O(E, d/i), 
where ^^©(E, dji) denotes the Hilbert space of holomor- 
phic functions which are square-integrable with respect 
to d/i on S. Due to this assumption U is an anti-unitary 
map and the holomorphic functions 



U|n) = = (T„(z), 71 = 0,1,... 

form an orthonormal basis in LP'0{Tj, dfi). 
One has the commutative diagram 



(4.42) 



U 



L^{R,da) L^O{J:,dp) 
where the anti-unitary map U o U^^ is given by 

(UoU-V)(2)= [ e~"'^lf{Lj)da{Lu), (4.44) 



(4.43) 



where ijj E (R, da) . Thus the Hamiltonian is given by 



(UoH/oU V)(z) EE i-^^'(z) 



(4.45) 



and is defined on the domain D{IJ o H/ o U 
e L20(E,d^) : £^ eL^O{J:,dfi)}. 
In terms of the Hilbert space L^O{T,, dpi) it is possible 
to find an explicit form of the creation operator a* , i.e. 
Hermitian conjug ate of the spectral annihilation opera- 
tor a defined by ( 4.31 ). W e wil l call a* th e spectral 
creation operator. Using (4.31) and (4.33), we obtain 



(Uoa*oU *)(z) EE z«'(z). 



(4.46) 



Thus we see that the domain of U o a* o U is given by 

^(Uoa* oU"^) 

= e L^O{T., dp) : e L^OiY., dfi)} . (4.47) 

Taking into the account the above considerations let 
us notice that the operator a* is described explicitly in 
the L^OCSjdiJ.)- realization and the operator a is ex- 
plicitly given in L'^(R,da)- realization. They satisfy the 
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canonical commutation relations 
[Ui,a] = [Uj,a*]^i 

with the Hamiltonian H/, giving 
[H/,a-a*] 



(4.48) 



(4.49) 



i.e. the operator a — a* is an integral of motion for the 
system under consideration. 



From a physical point of view (see (3.34, 3.35, 3.36|)) it 



is important to describe the time evolution in interaction 
picture of the system, i.e. e~'^^* 1-0(0)). To do this let us 
introduce the following notation for the matrix elements 
of e"*H^* 



Note that 



(t) (m|e-'«^*n)- / e-'"*P„(c^)P„(c^) da(w). 

(4.50) 



(4.51) 



where dnit) are given by (4.23) and they satisfy 



^n{t) = ^ (n|e-*H^*0) = P„(z -^)a(t). (4.52) 
The interaction evolution in the space J- is thus given by 

oo 

e-'-^'^m = ("^IV')^m,„(i)l"), (4.53) 

m,n— 

while in the space L^(M, dcr) the evolution is described 
by {i.i). In ^^©(I], d/x)-realization we have 



{\] oe^'^'^ o\] ^^^!{z) = ^{z + t). (4.54) 



As a consequence of (4.50) and (4.54) we obtain the re- 
lation 



'Sm,k{zi)dk^n{z2), 



(4.55) 



which for m = can be expressed in the form 

dn{zi+Z2) = 'Y'5k,7i{zi)dk{z2). (4.56) 



Moreover, putting m = n = in ( 4.55 ) we obtain the 
formulae (4.38) for the reproducing kernel. 



At the end of Section 3 it was shown (see( [3.34 ), (3.35) 
( |3.36 )) that the quantities 



(V(0)|e^"^'Fe-'"^*V(0)) 



(4.57) 



plays an important role if we consider the expectation 
values of the operator F on the time evolving state |V'(i)) 
(see (p8|)). 



The variety of the realizations of our model, namely, 
the L^(M, dcr) and L^OCS, dfi) repre senta tions allow 
us to give three equivalent formulae on (4.57) 



(V'(0)|e'"^*Fe-*"^*V(0)) (4.58) 

= ^ {mHK,nmn\Fk)dkAm\m) 



'm,'n,k,l 



e-'-^'^/jiuj) (U o F o U-i) {e"^'ip{uj)) da{uj) 



^(z + t) UoFoU (1'(x + i))d^(z,z), 



where = U|0(O)) and = XJ\ip{0)). In this way we 
have a very strong instrument for calculations of many 
physical characteristics of the system under considera- 
tion. 

In particular we have 



„iH/t 



ae 



-'Hit =a + t 



(4.59) 



and therefore 



0(O)) = (0(O)|a0(O))-t-i(V(O)|V(O)). 

(4.60) 



V. INTEGRABLE SYSTEMS RELATED TO 
CLASSICAL ORTHOGONAL POLYNOMIALS 



Here we shall investigate the classes of the physical 
systems with Hamiltonians of the form (|3.30D with the 
coefficients b{n) and h{n) given in table \ni[ The three 
classes of Hamiltonian operators are related to Hermite, 
Laguerre and Jacobi polynomials. We choose one mode 
case for simplicity and the circumstances which make the 
reduction not necessary. Then the Hamiltonians are ex- 
pressed in terms of usual creation and annihilation oper- 
ators in the following form: 



H 



Lag 



+ \/ a + a , 



h bo 2bi ^ bi 

— I ^ ^ 

ai Oi ai ai 



bi 



ai 



1 a* 



(5.1a) 
(5.1b) 
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Jac 



2a*a(a + b){^ + - 1) + 2{sL*a.y{a + b)- 2bfj. - 2av + /ii/(a + 6) + b^? + av^ 



+ (6-a)l 



+ (6- a)^ 



+ - 2 + 2a*a) + + 2a*a) 



(/^ + a*a)(i/ + a*a)(/i + + a*a — 1) 
{^jL + v + 2a*a - 1)(^ + V + 2a*a)2(^ + + 2a*a + 1) 

' (At + a*a + l)(;/ + a*a + l)(/i + ;/ + a*a) 

(/X + 1/ + 2a*a + 1)(^ + V + 2a*a + 2)2(^ + i/ + 2a*a + 3) 

I 



(5.1c) 



The ranges of the parameters v, bo, and oi are chosen 
such that the operators are well defined and are essen- 
tially self adjoint. Ini^(R, c?(t) (i.e. spectral) representa- 
tion the form ulae (3.3C ) lead to three-therm recurrence 
relation ( |4^ ) (see also ( A.lOj )). 



From Pearson equation (see (A_^) and table |) we ob- 
tain the expressions for measures: 



da^^'icj) = C'e'^*'^+^'' duj 



for cj e M, 



e ""i dio 



forc^e [-j^,^), 



(5.2a) 



(5.2b) 



(5.2c) 



for uj € (a, b) . 

In the holomorphic 
all Hamiltonians act 
formulae ( |4.45| )) but 
systems is hidden in 
diJ.{z,z) = fi{y)dxdy, ( 
of the domain S. The 



representation L'^0{'S, d^) 
as derivations: i-^, (see 
the difference between the 
the reproducing measures 
[z = X + iy), and the choice 
general case is described in 
Proposition IV. 3. Here we solve equation (4.36) for 
fi{y) in the special class, namely: continuous functions 



J 



except, possibly finite number of points in every compact 
subset. The discontinuity points are assumed to be of 
first kind. Let us summarize the results in the following: 



H) Hermite case: S = C and 



ai 
26n7r' 



(5.3.a) 



L) Laguerre case: 

i; = |z = .T + iyeC:2/<-^| and for /i > 1 



bQai 

2 



-2y 



bi 



2bQ 

e ""i ' 



(5.3.6) 



For fj, = 1 wc obtain an isomorphism of L^C'(E, d/i) with 
H^{D, dX) - the Hardy class of functions on the unit disc 
D C C with the measure dX supported on the circle dD = 
{e'"^ : (p e[0, 27r]} and given by 



dX 



1 



-dip. 



1 — sin ip 

J) Jacobi case: E = C and for 1.1 + u > 3 



(5.4) 



2 (6-a)i 



M^-(,) = } ^ r(M-l) 



C Tiv-l) ^ y' 



-(b+a)a(22^)^-2yy_^ 3_^[2(5_a)y] for y>0 

'^Wti^ 3_!i+.[-2{b- a)y] fory<0 



2 ' 2 2 



(5.3c) 



r 



where Wk,\{z) are confluent hypergeometric Whittaker's domain of parameter fj,, namely for fi > I, 
functions (for definition see [^|). This formula simpli- 
fies in the case ^ = v corresponding to the Gegenbauer 
polynomials. The following statement is true for a larger 



2 -2y 



C \b 



X e 



-(b+a)y_ 



1 



r(M-i)V^ 



K^.{{a-b)y){b.b) 
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with Ka{z) being the modified Bessel functions (for def- 
inition see jj]). 

For all three cases one can fi nd th e explicit form of 
matrix el emen ts of propagator ( 4.5C ). Because of the 
relations (4.52) an d (kl.51 ), we should display the charac- 
teristic functions (t4.9|) first: 



ai 



(5.6a) 



a^''9(z) = Cr(^)e 'i 



e-i'\ (5.6b) 



The symbols ct''"'^(z; /z, v) are introduced in order to sim- 
plify the next formulae. Using the Rodrigues formula 
(see (A. 8)) we obtain the explicit form of ct„(z): 



ci'"'{ib2zY d'^'^iz; fi + n,i^ + n). (5.7c) 



biz \ Tip + n) ^ 



.z + i-, 



(5.7a) 
(z), (5.7b) 



\n ^Jac 



d-''"'{z) = dj{z;fi,i^)^C 



r(/.)r(z.) 



r(M + ^^) 

X (6 - ar+^-'e'^''' iFii^U (a - fe)»z).(5.6c) 



After a simple but tedious calculation we find: 



\/\ ai/ ^ (m - fc)!(n - fc)!fc! 



m+n min{m,n} 



~2fc 



(5.8a) 



Lag Lag im+n^ + m)T{fl + v) ^Lag / \ /)t 



p ( ^I.+k-n. 



lai 



(5.8b) 



^Jac ( ^\ 



JacJac yn+n /m^ /«N / -, Nfc+i T (fl + m)r {^i + n) r{v + m)r{v + u) 

™ " ^ ^ii; r(/i + m - fc)r(/i + n - r(j^ + fc)r(i^ + 

X d''°-''{z; fi + m + n ~ k ~ 1,1^ + k + I) . 



(5.8c) 



The physical quantities which are of great importance 
are: the Hamiltonians H and H/, the creation a* and 
annihilation a operators, and the occupation number op- 
erator N = a*a. In our case the operators A, A* are 
also important. They can be interpreted as the cluster 
annihilation and cluster creation operators. Similarly, 
the operators a and a* which are related to the spec- 
tral coherent states map ( 4.25| ) are interested, too. Their 
physical meaning is partially explained by the commuta- 
tion relations (4.48). They are related to Hi, a and a* 



in the following way: 



V-aibo 



{bo + bi H^"9)a^°9 



= — z 6i a*a + ii/6^(a*a + /i) a 



(5.9a) 



(5.9b) 



J 



Jac\ ^Jac 



^ &(2/i — v + a*a — 1) + a{2v — /i + a*a — 1) 
^ + 1/ + 2a*a - 2 



(6 - a)(-/i - -I- a*a - 1)(^ -I- a*a)(i/ + a*a) / /i + + 2 a*a + 1 , „^ 
~^ rr-T-, — ■ tttt^ — ■ — rr: \ / ■ — a . (5.9c) 



|62|(Ai + i' + 2a*a- l)2(^ + iy + 2a*a)3 
I 



{^i + v + a*& - 1) 



In the spectral representation of H/ the operator a is by different formulae: 
jiven for all the cases by i-^ but the conjugates are given 



(5.10a) 
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for ^ > 1, 



aiuj + oq — 61 
biLu + bo 



d_ 

duj 



(fli + 262)^ + ao - 62(0 + ' 
62(0; — a){b — Lo) 



(5.10b) 



d 

dhJ 

(5.10c) 



for /X, V > 1. 

In the h olom orphic representation the operator a* is 
given by (4.46), i.e. as the operator of multiphcation 
by the argument z. The operators a^'^^ , a^""^ and a"^""^ 
are pseudodifFerential ones and we shall not express them 
explicitly here. 



The occupation number operators N defined by (3.25) 
take in the spectral representation the following form: 



^Lag 



fli J du! ai did'^ 
ao\ d / bi 



(5.11a) 



bo 



^ ai J duj ~^ \.ai ^ fli / duj^ ^ 



For the Jacobi case N-^""^ we are able to write down only 
the relation 



N 



Jac /j^Jac 



^ — 1/ — Vj = (uj — a){b — Lo) 



diO^ 



+ [{—y, — v)uj + ^b + va]—. (5.11c) 

duj 

In the holomorphic representation N can be expressed 
as: 



■^Lag _ 



60 



bi 





] z + z'^ 








1 dz 




bo \ 




I — z z 


Ol 


ai J 



ai 

,^b,.\ d 

1 H iz ] z— 

ai I dz 



(5.12a) 



(5.12b) 



Now, we will present the expectation values on the 
following states, interesting from the physical point of 
view: 



i) occupation number states |n), n € NU {0}, i.e. the 
eigenstates of N, N|rt) — n|n); 

ii) Gaussian coherent states C € C, i.e. the eigen- 
states of a, a|^) = CIC); 

iii) spectral coherent states |z), z £ E, i.e. the eigen- 
states of a , a\z) = z\z). 

Using the operators U and U one can realize these states 
i n spe ctral or holomorphic representation, too (see ([1.2|), 
(^lo|) , (|43^ ), (l^l)). 

Of course the Hamiltonian H/ does not depend on time 
t and its mean values are given by: 



(H/)„ 
(H/)c 



n=0 



hin) 

b{n+l) 
y/n+1 



(C + C) 



1 d ^ z — z 

2dy ^ ^ 2i 



(5.13) 

,(5.14) 
(5.15) 



The indices n, C, z are related to the occupation num- 
ber eigenstates, Gaussian coherent states, and spectral 
coherent states, re spe ctively. The function b(n) and h{n) 
are given in table III and a(z) is presented in (|5.6|). 

The mean values of the powers of the occupation num- 
ber operator are given as follows: 



(N'(i))„=^|i?„,fc(i)l'^', 
(N'(0)c=e 



k=0 



_,-lc|- 



\rm\n\ 



(5.16) 



fc'a™,fe(t)CT„,fe(t), (5.17) 



—^Y.\a^[^z + t)\^n\ (5.18) 
a\z — z) ^ — ' 



where / e N. It is interesting to rewrite the last formula 
for every polynomial class separately 



'(N^-^'Xt))' 



\z+ty 



. aij 



iFi 



/-I 



.2 



(5.19a) 



(N^"5(t))' 



9^ z — z 

^61 2i 



\z + 



\z + t 



bi I 









f /j+1,2... ,2 


Z 4- t 


1 




z+t-i^ 

Ol 





(5.19b) 



\b^{z + i)|"" \a"'\z + i; M + n, J/ + n)\^n' 



J Jac 



(5.19c) 
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We give now the formulae for the correlation functions: 



^ (jn + r)!(TO + s)! 



to! 



(5.20) 



/c.m.n— 



(5.21) 



^ OO 

(a"(i)a"(i))^ ^ — -Y,dn+r{z + t)dr.+s{z + t) 

^ ' n— 



^(n + r)!(n + s)! 



(5.22) 



Replacing the creation and annihilation operators a* , a by the cluster creation and the cluster annihilation operators 
A* and A we obtain the functions which by analogy will be called the cluster correlation functions: 



(A"(t)A^(0)„ = Kk + r)...b{k + mk + s)...bik + l)an,k+rit)^k+sAt), 
k,i=a 



(5.23) 



1^,771,1 — 



'mill 



b{k + r)...b{k + l)b{k + S)...b{k + l)^77^,k+r{t)^k+s,l{t), 



(5.24) 



(A"(t)A'(t)), = ^^^^ £ bik + r)...b{k + l)b{k + s)...b{k + l)dk+r{z + t)dk+s{z + t) 



k=0 



(5.25) 



The time evolution of a is given by ( [4.59|) . This allows 
us to express the time dependence of (a (t))^, (where 
\ip) is an arbitrary state and Z € N), in terms of the mean 
values of some powers a = a(0) acting on the state \tp): 



fc=0 



(5.26) 



As a consequence we conclude that the dispersion 
(Aa(t))^ = Hm^]^) of the oper- 

ator a{t) in an arbitrary state \ip) does not depend on 
time 



the following strong- field Hamiltonians H^: 



Hi 



H 



Jac 



(a + a ), 

,^1 



Hf-^f = -2— a a 

ai 

bi 



ai 
a + b 
2 

b~a / 1 



(5.30a) 



(Va*a + la+\/a*a + 2a*), (5.30b) 



4 V\/a*a + 1 Va*a + 2 



a* .(5.30c) 



(5.27) 



The following expectations take especially simple form: 



(5.28) 
(5.29) 



These Hamiltonians belong to the respective families 
given by ( |5.l| ). They a re ob tained in the Hermite case 
by putting ao = in ( 5.1a ), in the Laguerre case by 



putting fJ. = 1 and bo = — ^ in (5.1b), and in the Jacobi 



case by putting /i = = | in (5.1c) 



Let us recall the definition of the phase operator [E3| 



Let us now see what will happen when the intensity of 
electromagnetic field is sufficiently large (the light of a 
strong laser). This cor respo nds to the limit of the large 
n in the Hamiltonian (3.30) (see table too). We get 



exp := (a*a + 1) ^a., 
exp {-i(f>) := a*(a*a + 1)~^, 

cos {(j)) :— i (exp {i(j>) + exp {—i<j>) 



(5.31) 
(5.32) 

(5.33) 
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We can now rewrite ( 5.301:1 ) and (5.30c): 

bi 



H 



Jac 



ai 

a + b b — a 



2a*a + 2a*a cos (0) + exp (i(/)) j ,(5.34) 
cos(^). (5.35) 



So, in the Jacobi case in the strong-field limit, the Hamil- 
tonian tends, up to a constant, to the cosine of the phase 
operator. This subcase does not depend on the choice of 
the ranges of the parameters /i, v. 



VI. A PHYSICAL REMARKS 



Parametric modulator 



In order to pre sent some physical interpretations of 
the Hamiltonian (2.25) with H/ given by ( ^.37 ), let us 
rewrite it in the following form 



M 



H/ = ^ (^ja-ja^ + h{a.Qao, ... , a^^ aAf) 



it 2 i^j^jS 



'g(a*ao, . . . , a^^aj,/)a[f . . . a^f/ 



it 2 li'ja's 



'5(a*ao, . . . ,aX^aM)ao° . . .a'" 



M 



(6.1) 



The first term, which is linear in photon number opera- 
tors describes the free field. The second term, which is an 
arbitrary function of these operators may be treated as 
a generalization of the Kerr medium description, where 

H/ = ^ ^(a*a)^ — a*a^ , where x is proportional to 

the third-order nonlinear susceptibility, [^. The terms 
of the type /i(aQao, . . . ,a^^aM), after the appropriate 
choice of the function h, play an important role in the 
theory of the nondcmolition measurmcnt ||l^ and 

in the description of many other phenomena e.g. the 
optical bistabihty effect Q. 

The last term in ( ^.l| ) one can interpreted as a gen- 
eral form of the parametric modulator Hamiltonian. To 
motivate this interpretation let us recall the form of the 
Hamiltonian of nondegenerate parametric amplifier [G7[, 



H = o^oaQao -t- ujialsLi 
+tg (e2-*aoai 



2iLjt 



aoai 



(6.2) 



This Hamiltonian describes the case when the classical 
pump mode at frequency 2lo interacts in a nonlinear op- 
tical medium with two modes at frequency loq and wi, 
such that ujQ -\- LOi = 2w. If the system starts in an ini- 
tial Gaussian 2-photon coherent state |CoCi)i the mean 
photon number in 0-mode after time t is 

(a^(t)ao(t)) = |Cocoshgt-|-Crsinh.gtp + sinh^gi,(6.3) 

hence this mode is amplified. The next example is the 
Hamiltonian for the frequency up-converter (|27]|) 



H = ti>oaoao + ^^lal&i + k (e^'^'ajaj^ 



where uj = loi — ujq. 



(6.4) 



It is easy to compare (6.1) with (6.2) and with (6.4) 
and conclude that our Hamiltonian is a natural gen- 
eralization of that describing parametric am plification. 
In order to understand that in general (6.1) describes 
not only amplific atio n but also modulation let us no- 
tice that due to (3.3) we can express the mean values 
(a* {t)aj (t)), j = 0, . . . , M in terms of the mean values of 
the operators Ao(t), .. . ,AM{t). But Ai(t),... ,AM{t) 
are the integrals of the motion, so if our system starts at 
the initial state from the reduced subspace T (see ( ^.26[ )) 
we obtain 



{a*{t)aJt))=l,{Ao{t))+/3„ 



(6.5) 



where the constant /3j are uniquely determined by 
Ai,... ,Xm and the matrix a. This means that the 
mean photon number in each mode is a linear function 
of (Ao(t)) or, in other words, the strength of the light 
in each mode is modulated by the function {Ao{t)). The 
modulation of the j-th mode depends on the exponent Ij . 
The shape of the function {Ao{t)) depends on t he c hoice 
of the coupling function ^(aQao, . . . , a^^aj\/) in (6.1) and 
the initial state of the system. 

As an example of the modulation function {Ao{t)) let 
us consider the situation when after reduction, we ob- 
tain the case corresponding to Laguerre polynomials and 
t he ini tial state is the spectral coherent state \z). From 
( ^.19b| ) we obtain 



(Ao{t)),^E\z + t\^ +F, 



(6.6) 



where the real constants E,F depend on /i, ai, bi 
and Ao,;. In this example the modulation function is of 
parabolic shape. This means, that in some interval of 
time we have the amplification and dumping of the light 
signal in others. 
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B. Generalized squeezed states 



The special cases of the interaction evolution operators 

-iHrt 



are the unitary displacement operators ||1 
D(C)=exp(Ca*-Ca), C e C, 
the unitary squeeze operators |^ 

S(z) = exp(za2 -za*2), zeC 
and the unitary two-mode squeeze operators 
T(e) = exp(eaoai - ^a*a.l), ^ e C. 



(6.7) 



(6.8) 



(6.9) 

This means that the Glauber coherent states and the 
squeezed states are special cases of the spectral coherent 
states defined in Section 4. In such a way, the two con- 
cepts of the notion of the coherent states meet each other 
in our framework. The first one, presented in |24[ , 
1^, is related to the minimalization of the suitable 
uncertainly relations. The second one, presented in [ pO[ , 
PH is based on the symplectic embedding of the classical 
phase space of the system into the quantum phase space 
(equipped with the Fubbini-Study symplectic form). 



and Q satisfies the Pearson equation with appropriately 
chosen polynomials {A{uj), B{uj)). Namely if deg B{lu) — 
(i.e. &2 = &i = 0) then we obtain the Hermite polyno- 
mials; if deg B{uj) — 1 (i.e. 62 = 0, 61 7^ 0), we obtain 
the Laguerre polynomials and if deg B{lo) = 2 (i.e. 62 7^ 
0), we obtain the Jacobi pol ynomials. In the last case 
the boundary conditions (A.4) hold if and only if a and h 
are roots of B{uj). For solution of the Pearson equation 
in these cases see table |. Additional conditions enforced 



on A{uj) by (A.4) are presented in this table too. 

By straightforward calculation one can prove that the 

-^P„(w) > , /c e N, is orthog- 
onal in the space L2(M,dCT(''')), where 



da'^'^^Lo) := B''[Lo)da{Lo). 



(A.6) 



The weight function g'^^^ satisfies the Pearson equation on 
the interval (a, 5) associated with {A^^\ijj), B{uj))^ where 



A'^'^^uj) ■.^A{uj) + k 



dB(uj) 
d(jj 



(A.7) 
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Appendix 

Here we present some facts from the theory of classical 
polynomials. 

Let us consider a pair of real polynomials {A{lu), B{lo)) 
of degree not greater than one and two, respectively 



A{u:) 



aiUJ + oq. 



-8(0;) := 1)21^^ + biuj + bo 



bi e 



(A.l) 
(A.2) 



The Pearson equation associated with {A{lu),B{uj)) on 
the interval (a, 6) C M (—00 < a < b < +00) is the 
differential equation for the weight function g: 



— [gB) = gA 
duj 

with the boundary conditions 
g{a)B{a) === = g{b)B{b). 



(A.3) 



(A.4) 



Each family of classical orthogonal polynomials {Pn] can 
be obtained by the Gram-Schmidt orthogonalization of 
the basis {w"}J^q in i^(R, dcr), where 



Lo < a 

da{uj) := ^ dto a < u <b 
w > 6 



(A.5) 



Proposition A.l For a given Pearson data i.e. a pair 
{A{lj)^ B{uj)) on {a,b) C R the following statements are 
equivalent: 

A. {P„([j)}5^Q form an orthonormal system in 
L^{R,da). 

B. The polynomials are given by Rodrigues ' formula 

1 r/" 

Pn(a;)=c„^.— (e(c.)i3"(u;)), (A.8) 

where c„ is the normalising constant. ( see table ^ 

C. The polynomials {-Pn('^)}^o satisfy the differential 
equation 



"^^^^d^ ^ ^^^^IhS^ ' ^"^"^^ " A„P„(w), (A.9) 



where A„ — ain + b2n(n — 1). 

D. The polynomials {P„(cL;)}^g are related by the tree- 
term recurrence formula (for h{n) and b{n) see table 

m 

uPnitu) = h{n)Pn{uj) -f 6(n)P„_i(w) + b{n + l)P„+i{uj) 

(A.IO) 

with the initial condition 



Po{uj) = const 



da{ijj) 



(A.ll) 
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TABLE I: 

Pearson data Weight function Additionally conitions 

A{u!) B{uj) (0,6) 

= oiw + 00 -B"'=''(a;) = 60 (-00, 00) q"''''{uj) = Ce^'"'+^^' C > 0, < 

yl^'^'^H =aiw + ao = 6iw + 60 (-|^,oo) ^'^"''H = C(u; + ^)'^-ie^" C > 0, fJ- < 0, 

H := > 

C > 0, o < 6, 62 > 0, 

^•^"=H = aiw + ao B''"=H = 62(w-o)(6-a;) (0,6) ^-^"=(0;) = C(a; - "^(6 - o;)""^ M := fffs^ > 0, 

:= ^2i±iiii > 
62(01—6) 

TABLE II: 

Cn — in Rodrigues' formula 

c«- = (Cn!(-ai&o)V-'^f')"^ 

= i^Cnl (-aibir (-^) r(/x + n)e ^ j 
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b(n) 








(n) 


= ^J-> 






^ ai 






= -^v"i 


n+fi - 1) 












n(/x+n — l)(iv-|-ri- 


-l)(lJ. + J^ + n-2) 




-l) + 2n^ (a + b)-2b^i-2av + ^iv{a+b)+bfi^ +av'^ 








n,-2)^(;i+i^ + 2n-l) 




{fi+v + 2n — 2^)(ii+iy+2n) 
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